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ABSTRACT 

Context. The interaction of solar oscillations with near surface convection is poorly understood. These interactions are likely the 
cause of several problems in helio- and astero-seismology, including the so-called surface effect and apparently unphysical travel time 
shifts as a function of center to limb distance. There is thus a clear need for further theoretical understanding and observational tests. 
Aims. The aim is to determine how the observed modes are affected by the convection. 

Methods. 1 use HMI velocity and intensity images to construct k-ai diagrams showing how the oscillation amplitude and phase 
depend on the local granulation intensity. 

Results. There is a clear and significant dependence of the observed properties of the oscillations on the local convection state. 
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1. Introduction 


The effect of convection on wave propagation is a complex is¬ 
sue. With no clear separation of scales in time or space, simple 
approaches do not work well. This is particularly unfortunate as 
there are observed effects for the Sun and other stars which are 
suspected of being due to such effects. While attempts at more 
comp l ete treatments have been made (e.g. St x&^hugzhdal 
I 2 OO 4 I: iHanasoge et al.llT013l : iBhattacharva et ^ 2015h . the de- 
tails are still far from understood. 

One example is the so-c alled surface effect, which has been 
known for a long time (e.g. IChristensen-Dalsgaard et aklllQSSll 
and which manifests itself as a difference between the theoretical 
and observed mode frequencies and has properties su ggesting a 
near su rfac e origin. For r e cent d iscussions of this see lPiau et alJ 
Jml and iBall & Gi^(l2^ . 


Another troubling effect is a time shift observed in time- 
distance helioseismology with an apparently unphysical depen¬ 
dence on the observable an d a strong center-to-limb dependence 
(se e. e.g. |Zhao et aT| 2012L An attempt to explain this was made 
bv iBaldner & Sch~ (l2012l) . who noticed that eigenfunctions in 
hydrodynamic simulations exhibit strong phase variations with 
depth and were able to roughly explain this, as well as the center- 
to-limb effects, using a crude theoretical model incorporating 
some of the properties of the convection. 

Here I address the latter of these two issues by investigating 
how the appearance of the oscillations depends on whether one 
observes them in the middle of granules or in intergranular lanes. 


2. Method 

While the interaction of waves with convection is complex, it is 
probably reasonable to try a simple model and start by consid¬ 
ering a horizontally uniform and initially vertically propagating 
wave in a Cartesian box. The two largest contributors to the dis¬ 
tortion of the observed wave are likely the local thermodynamic 
state (i.e. sound speed, as given by the temperature and thus in¬ 
tensity) and the vertical flow, which are highly correlated. Here 


I use the intensity, as it is easily measured and is less affected by 
the modes than the Doppler velocity. It is probably also reason¬ 
able to linearize the variations, that is to assume that the observed 
velocity varies linearly with intensity. Of course, the effect is not 
expected to be a real multiplicative factor, nor independent of 
frequency, so is best described in the Fourier domain; 


V(io,I') = Vo(co)(l + a(tj)I') 


( 1 ) 


where ~ indicates Fourier transform, V is the line-of-sight ve¬ 
locity, Vo the unperturbed velocity, /' is the intensity minus the 
mean (over time and space), ai the frequency, and a a complex 
factor to be determined. 

To measure the effect I us e data cubes from t he Helioseismic 
and Magnetic Imager (HMI: l^hou et al]l2012h of Doppler ve- 
locity (V) and compute d continuum i ntensity (/) , track ed by 
iLangfellner et al.l (1201 4l) at the rate of I Snodgrass! (1 19841) eval¬ 
uated at the center of each cube. These have 512 by 512 pixels 
of size 0.0005 radians by 0.0005 radians (» 348 km by 348 km) 
and 1920 points spaced by 45 s in time, for a total cube size of 
« 180 Mm by 180 Mm by 24 hours. For the present study cubes 
centered on latitudes 0°, +20°, +40° and ±60° crossing the cen¬ 
tral meridian near the center of the time interval are used. A few 
missing frames each day are filled by interpolating the adjacent 
images linearly. Also, for ease of interpretation, /' = /// - 1, 
where 7 is the mean (over time and space) of /, is used. 

Each frame of the cubes is subdivided into 4x4 pixel areas 
and two quantities are calculated; V, which is the average of V, 
and S which is the slope from a linear fit of V versus /. Rather 
than oversampling, the resulting images have 1/4 as many pixels 
in each direction. 

In reality a might depend on the horizontal wavenumber of 
the wave, so instead of simply averaging horizontally, the result¬ 
ing cubes are circularly apodized between 0.85 and 0.95 of their 
half width with a raised cosine and Fourier transforms V and S 
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are calculated from V and S. From these, a is obtained from the 
ratio of the cross spectra and power spectra as: 

{S{k^,ky,u,i)V{k^,ky,u),iy) 

a{k,oj) = -, (2) 

{y{k^, ky, oj, i)V(kx. ky, u, i)*) 

where () denotes averaging over azimuth 0 and time (indexed 
by the cube number /)[J, k^ = kcoscp, ky - ksin^, k - 1 /Rq 
is the (horizontal) wavenumber, I the spherical harmonic degree, 
* indicates complex conjugation, and Rq the solar radius. For the 
results shown here the time averaging was over 30 one day cubes 
starting on 2010 June 1. 

When determining V and S it is effectively assumed that the 
horizontal wavelength is long compared to 4 pixels, which may 
lead to some degradation of the results near the spatial Nyquist 
frequency of the V and S maps. However, far from the Nyquist 
frequency this should not be a problem. Similarly, the use of V 
as an estimate of Vq should be a good approximation as the term 
a{<jj)I' in Eq.[T]is small. 


3. Results 

Fig. [2 shows that there is a strong signal in both the ampli¬ 
tude and phase of a and that the modes behave quite differently 
from the background granulation signal. To get a better picture. 
Figure |2] shows the amplitude and phase as a function of fre¬ 
quency along the ridges. It is clear that both change dramatically 
as a function of frequency. This is, perhaps, not surprising since 
the frequen cy determines the near s urface phase behavior of the 
modes (see lBaldner & Schoull2012ll . 

However, there is still significant scatter between the differ¬ 
ent radial orders, especially for the phase. Some of this is un¬ 
doubtedly due to the varying oscillation/convection power ratio, 
which in itself makes the very low and high frequency results un¬ 
reliable. Another thing to consider is that the oscillations are vis¬ 
ible in both V and I. Indeed, studying V and I together provides 
information on the mode exc i tation and mode physics, as de - 
scribed in, e.g.. lSeverino et all(l2008h and lSeverino et al.l(l2013l) . 
This also means that some of the effects seen may be due to the 
visibility of the modes in /, which is used to determine S. To 
correct for this, I' is low pass filtered with a cutoff of 2 mHz, 
to eliminate most of the mode signal and the analysis repeated. 
As Fig. [2 shows the scatter of the phase is indeed significantly 
reduced. 

The amplitudes are quite interesting. The change with inten¬ 
sity is close to zero at the lowest reliable frequencies (i.e. » 2.5 
mHz) but becomes larger than unity above about 3.5 mHz, in 
other words a unit change in background intensity (correspond¬ 
ing, roughly to 100 km height change) causes of order unity 
change in the amplitude in the observed velocity. 

It is also worth nothing that the phases are generally not close 
to 0 or 180 degrees. In other words the phases at different inten¬ 
sities are not generally in phase. 

These results are , perha ps, not surprising given the results 
of iBaldner & Sch^ (l2012h . who predicted large phase varia- 
tio ns with height. A quant itative comparison is quite difficult, 
as IBaldner & Schoul (12012h did not address the same problem 
and since their final result depends on a delicate competition be¬ 
tween positive and negative contributions from different heights 
and details of the radiative transfer. But noting that a unit change 

* The averaging over (p and time is done by first interpolating each of 
the cross- and power-spectra to a grid in k and 0, averaging over 0 and 
then averaging over time. 



I 


Fig. 1. From top to bottom, (a) shows the velocity power on a 
logarithmic scale with a range of 10^ (blue-green-red-yellow), 

(b) the magnitude of a on a logarithmic scale from 0.1 to 10 and 

(c) the phase of a. The color coding for the phase is shown in the 
inset. Green (0°) corresponds to the amplitude increasing with 
I, black (180°) to decreasing, blue (90°) to the waves at higher 
intensity leading those at the reference intensity and red (270°) 
to trailing. Rough mode frequencies (courtesy of R. Bogart) for 
radial order 0 < « < 10 are shown by dots. 


in the relative intensity corresponds to roughly 100 km height 
difference and assuming that there are no other effects, they 
saw phase changes of order tens of seconds (of order a radian), 
so it is not surprising that we observe changes of order unity. 
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Fig. 2. Cuts through panels b and c of Fig. [1] along the ridges 
for «=0 through 4. Dotted black line shows the results for the 
noise evaluated half way between the ii - I and n — 2 ridges. 
For the phase convention see the caption to Fig.[T] Note that the 
large oscillations in the phase just below 2 mHz are caused by the 
very low amplitude in that region and thus large noise sensitivity. 
This, in turn, also results in problems trying to unwrap the phase, 
which might otherwise have improved the presentation. 


However, more detailed predictions would require more detailed 
theoretical work. For at tempts at this the reader is referred to 
IStix & Zhugzhdal (l2004h and references therein. Similarly, sub¬ 
stantial insight might be obtained from numerical simulations. 

As for the robustness of the results, it may be noted that us¬ 
ing the low pass filtered Doppler velocity instead of the intensity 
leads to very similar results, except for a sign change (veloc¬ 
ity and intensity are strongly anti-correlated in the granulation). 
Also, assuming that the variation is linear with intensity and that 
the velocity measurement process is linear, it follows that the 
slope is insensitive to blurring and thus to the instrument PSF 

As a likely cause of the phase shifts is the height dependence 
of the phase of the oscillations combined with the varying height 
of observation, one might expect some change with viewing an¬ 
gle, as that also changes the effective height of observation. As 
seen from Fig.|4]the phases are almost unaffected up to 40° from 
disk center. At 60° the results are dramatically changed, but are 
different for -1-60° and -60°, so this is likely not a robust result, 
most likely due to the substantial foreshortening. The amplitudes 
show more consistent changes. To what extent those changes are 
due to changes in the S/N has not been investigated. 



Fig. 3. As Fig.|2] but using intensities passed through a low pass 
filter with a cutoff frequency of 2 mHz. 


One might also expect the intensity signal of the oscillation 
to change with position within the granulation pattern, which as 
Fig. |5] shows is indeed the case. Here the raw intensities were 
regressed with the filtered ones, in the same way as was done for 
the velocities. Given the lower S/N for intensity, these results are 
probably more affected by noise, but it is still clear that there is 
a strong signal and while the results for the modes and the back¬ 
ground appear similar around 3 mHz, the S/N is actually quite 
significant there. It is interesting that the results are substantially 
more n dependent, probably reflecting differences in the mode 
physics (e.g. f modes do not compress the matter significantly, 
while p modes do). It is also interesting that the frequency depen¬ 
dence is different from that of the velocity, presumably reflecting 
the different heights of formation. It has not been investigated 
how other variables, such as line depth or width, vary. 

Finally it may be noted that previous studies of waves over 
granul es and intergranu l ar lan es h ave been made (see , for ex¬ 
ample, l^omento^et^ (1200 ih and iKostvk et al.l (l2006h and ref¬ 
erences therein). Those studies were mostly concerned with at¬ 
tempts at detecting the excitation of modes or interactions with 
magnetic fields. Also, they did not investigate the dependence on 
k and a> and so were not able to resolve the modes clearly, which 
makes comparisons with the present work difficult. 
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Fig. 4. Amplitudes and phase shifts for « = 2 at various lati¬ 
tudes. Note that the average BO angle was only about 1°, so the 
viewing angle is effectively identical to the latitude. 


4. Conclusion 

Using the method described it is shown that the amplitudes and 
phases of the oscillations depend strongly on position in the 
granulation. Of course, this represents a correlation and does not 
imply that the intensity is the proximate cause. The velocity and 
intensity are highly correlated and no attempt was made to dis¬ 
tinguish between their effects. No physical model of the effect 
has been m ade, but the large phase ch ange with geometric height 
reported bv iBaldner & Sch^(l2012h combined with the change 
of observing height with intensity in the granulation is a promis¬ 
ing candidate. Also, instrumental effects can not be completely 
excluded, especially for the amplitude effect, as there is some 
change in Doppler sensitivity with the thermodynamic proper¬ 
ties where the line is formed. 

A promising way to investigate the origin of the effects re¬ 
ported here is to analyze the results of large scale numerical sim¬ 
ulations of near surface conv ection. Indeed, studies similar to 
those of ISeverino et al.l (l2013l) . but for the effects studied here, 
may prove very useful. 

Observations with higher spatial resolution may also help 
shed more light on the relevant processes. In particular it may 
be possible to look for non-linearities and to disentangle the in¬ 
tensity and velocity effects. Similarly, higher spectral resolution 
may allow for the height dependence to be determined. It may 
also be possible to correlate against other variables, such as the 


Fig. 5. As Fig. |3]but for the intensity. 

magnetic field strength, alone or in combination with the inten¬ 
sity, thereby gaining further insight about e.g. magnetic field ef¬ 
fects, a subject of significant current in terest. _ 

Note that, similar to the analysis bv iNagashima et all (l2014ll . 
the analysis here only relies on existing HMI observations and 
does not require new observations using higher spatial resolu¬ 
tion or other spectral lines. It is thus possible to study and ex¬ 
ploit the effects described here using the existing 5 years of HMI 
data. An interesting possibility is that the effects described here 
may be exploited to further separate granulation from oscilla¬ 
tions, thereby possibly leading to an improved S/N ratio. 
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